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I investigate entanglement entropy of one dimension (ID) cluster states to discuss the delocal¬ 
ization of clusters in linear-chain 3a- and 4a-cluster states. In analysis of entanglement entropy of 
ID Tohsaki-Horiuchi-Schuck-Ropke (THSR) and Brink-Bloch cluster wave functions, I show clear 
differences in the entanglement entropy between localized cluster wave functions and delocalized 
cluster wave functions. In order to clarify spatial regions where the entanglement entropy is gener¬ 
ated by the delocalization of clusters, I analyze the spatial distribution of entanglement entropy. In 
the linear-chain 3a cluster state, the delocalization occurs dominantly in a low-density tail region 
while it is relatively suppressed in an inner region because of Pauli blocking effect between clusters. 

In the linear-chain 4a state having a larger system size than the linear-chain 3a state, the delo¬ 
calization occurs in the whole system. The entanglement entropy is found to be a measure of the 
delocalization of clusters in the ID cluster systems. 


I. INTRODUCTION 


A variety of cluster states have been known in light nuclei, such as 2a states in 8 Be, 16 O + a states in 20 Ne, and 
3a in 12 C (for example, see Ref. [l] and references therein). In this decade, a new concept of cluster states has been 
proposed to understand cluster motion in these states @t 3|- That is a dilute cluster gas state where clusters are not 
spatially localized in certain positions but they are rather freely moving like a gas. To describe the non-localized 
(delocalized) cluster states, a new type of cluster wave function, so-called “Tohsaki-Horiuchi-Schuck-Ropke” (THSR) 
wave function, has been introduced. The THSR wave function is essentially based on a clusters in a common Gaussian 
orbit having a range of the system size. It has been shown that the 2a state for 8 Be(0) 1 ’) and the 3a state for 12 C(0^~) 
can be well described by the THSR wave functions with a large Gaussian range compared with the cluster size, and 
therefore, these states are interpreted as gas-like cluster states of 2a and 3a. Recently, Zhou et al. have extended 
the THSR wave function to apply it to 20 Ne system, and shown that the THSR wave function can also describe the 
16 0+a states in 20 Ne [fjj, 0]. 

More recently, Suhara et al. have proposed that this concept of the a-cluster gas is applicable also to one dimension 
(ID) cluster motion in linear-chain na structures ||. Existence of the linear-chain na states has been a long standing 
problem. In the early stage, Morinaga proposed a linear-chain 3a structure in 12 C(0^") (sj, but this assignment was 
excluded from the experimental data of its a decay width [Tol |. The possibility of the linear-chain 3a structure in a 
hi gher 0 + state is now on discussion. Negative results for the linear-chain state were obtained by 3a-cluster models 
0,11|, whereas appearance of a chain-like 3a state with an open triangle configuration was predicted by microscopic 
approaches with no_cluster assumption and by an approach with the cluster breaking effect |l2H15j and also by an ab 
initio calculation fl6| . Also linear-chain 4a structures in 16 O have been attracting a great interest, and searching for 
linear-chain states have been performed in experimental and theoretical works in this line fl7l - [22j . 

The conventional picture for the linear-chain na structures is spatially localized a clusters arranged in ID with 
certain intervals. Such a localized na state is described by a single Brink-Bloch (BB) wave function (23). The 
essential result shown by Suhara et al. is that the wave functions of the linear-chain 3a and 4a states are described 
by superposing a large number of BB wave functions, and surprisingly, they have large overlaps with 1D-THSR wave 
functions Q. This result indicates that the delocalization occurs in ID a-cluster motion and the linear-chain states 
can be regarded as ID a cluster gases of delocalized clusters rather than localized cluster states. However, as shown 
in Ref. [8j, the density distributions of the linear-chain states of 3a and 4a have three- and four-peak structures, 
respectively, which indicates partial localization of a clusters because of the Pauli blocking (repulsion) between a 
clusters. Even though the model in Ref. Q was restricted in ID configurations and the stability of the linear-chain 
states against bending motion and a decays have yet to be investigated for a conclusive answer to existence of the 
linear-chain states in realistic nuclear systems, their work provides the new picture of ID cluster states, which is 
important to understand cluster phenomena in nuclear many-body systems, and also academically interesting. 

In light nuclei, clusters are often formed by spatial correlations of constituent nucleons. Even in a mean-field 
picture, spatially correlated nucleons in a cluster can be described by a product of localized single-particle wave 
functions. However, such a localized cluster usually has much kinetic energy for the localization of center of mass 
motion (c.m.m.) of the cluster. If there is no potential nor Pauli blocking effect between clusters, it is naively expected 
that delocalization of the c.m.m. of clusters occurs to release the kinetic energy. An ideal case of the delocalization 
limit is a zero-momentum cluster gas state. The delocalization of clusters in realistic cluster states may depend on 
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competition between kinetic energy gain and energy loss by other effects such as potential and Pauli blocking between 
clusters. 

To distinguish between localized and delocalized cluster structures in microscopic wave functions, one should care¬ 
fully consider the antisymmetrization of nucleons between clusters, which strongly affects inter-cluster motion at a 
small distance. The antisymmetrization effect suppresses an amplitude of the inter-cluster wave function at a small 
inter-cluster distance and works as the Pauli repulsion (blocking) between clusters, whereas the effect vanishes at a 
large distance. If the system size is comparable to the cluster size, clusters can not move freely because of the Pauli 
blocking effect. Therefore, the delocalization of clusters likely occurs not in high-density regions but in low-density 
regions. Indeed, for two-body cluster states of 2 a and 16 0 + a in 8 Be and 20 Ne, the author has shown that the delo¬ 
calization of a clusters occurs at a long tail part of the inter-cluster motion [24j. Also in the result of the linear-chain 
3a and 4a states shown by Suhara et al ., the peak structure of density distributions shows the partial localization 
of a clusters in an inner region as mentioned previously. It may suggest the possibility that the delocalization does 
not occur in the whole region of the system. I should stress here that the THSR model wave functions, which can 
successfully describe cluster gas states, have characters of localization and/or delocalization of clusters depending 
on the system size relative to the cluster size. In the case that the system size is large enough compared with the 
cluster size, the THSR wave function actually describe a dilute cluster gas. However, if the system size is as small 
as the cluster size, the THSR wave function become equivalent to a localized cluster wave function because of the 
antisymmetrization. A general question is how one can understand the delocalization of clusters in an intermediate 
case between both limits of localization and delocalization. To clarify the region where the delocalization occurs, one 
may encounter a difficulty from the antisymmetrization, i.e., Pauli blocking effect. The Pauli blocking effect often 
makes it difficult to tell the difference between localization and delocalization of clusters in a relatively high-density 
region, where the delocalization is usually suppressed. Note that one can not obtain a definite answer from a wave 
function without the antisymmetrization, which may contain unpliysical forbidden states. 

My aim is to analyze the delocalization of clusters in microscopic wave functions with an approach free from 
the antisymmetrization effect. For this aim I propose a method of analysis using entanglement entropy defined by 
the one-body density matrix. The entanglement entropy has been introduced by Bennett et al. in 1996 25j], and 
widely used in various fields such as condensed matter and quantum field theory (see for example Refs. [26H28| and 
references therein). The entanglement entropy indicates how particles are entangled with other particles, and can 
be an indicator to measure many-body correlations as applied in in various fields such as condensed matter physics 
and quantum physics. For a system of independent Fermions, the wave function is given by a Slater determinant 
of single-particle wave functions, and the entanglement entropy completely vanishes. The entanglement entropy is 
generated by many-body correlations beyond a Slater determinant. In the case of cluster states, the entanglement 
entropy is zero in a single BB wave function for the localized cluster limit, and it is generated by the delocalization 
of clusters. 

In this paper, I analyze the entanglement entropy in ID cluster states for na and 16 0 + a systems. Based on the 
analysis of the entanglement entropy, I investigate the delocalization of clusters in the linear-chain 3a and 4a states 
predicted by Suhara et al. as well as the 2a state in 8 Be(0^) and the 16 0 + a states in 20 Ne(0/") and 20 Ne(lj“). 

This paper is organized as follows. I describe the formulation in Section [TTJ Section Hill discusses the entanglement 
entropy of ID cluster states. The paper concludes with a summary in section HV1 


II. FORMULATION 

A. Wave functions for linear-chain na-cluster structure 

I consider the linear-chain na-cluster structures aligned to the z axis (n is the number of a clusters). For simplicity, 
the angular momentum projection is not taken into account, and only ID configurations of n a clusters in intrinsic 
wave functions are considered in the present analysis. It means that the (de)localization of a clusters are defined for 
a-cluster motion along the z axis. In this section, I first describe a general form of a-cluster wave functions for the 
linear-chain structures. More details of practical model wave functions used in the present paper are described in the 
latter part of this section. 
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1. BB wave function 

I use the BB wave function [23| for a localized na-cluster wave function, 


^(R u ...,R n ) 


(1) 
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is the four-nucleon wave function of the ith a cluster expressed by the (Os ) 4 harmonic oscillator (ho) shell-model 
configuration localized around the spatial position R,. Note that a single BB wave function for an na system is 
written by a Slater determinant of single-particle wave functions. For a linear-chain structure aligned to the z axis, 
the position parameter i?,; is set to be Ri = (0, 0, Ri), and the ID BB wave function is expressed as 3 >bb(-Ri,..., R n )- 
The parameter b for the a-cluster size is chosen to be b = 1.376 fm same as in Ref. Jj. General wave functions for 
ID na systems can be written by linear combination of BB wave functions $ 33 (^ 1 ,..., R n ). 


2. 1D-THSR wave function 

As shown in Ref. the linear-chain 3a- and 4a-cluster states in 12 C and 16 0 systems are well described by the 
1D-THSR wave functions proposed by Suhara et al.. The 1D-THSR wave function is given by linear combination of 
BB wave functions with a Gaussian weight and it is in principle written as 



where Xi is the center of mass (c.m.) coordinate of the iih a cluster and (f>(ai ) is the intrinsic wave function of the a 
cluster. If the antisymmetrization is ignored, the 4 >( ! ' 3 -thsr(^) expresses the na state where all a clusters are confined 
in the x and y directions in the size b/y /2 while they move in the z direction in the Gaussian orbit specified by the size 
parameter /3, which corresponds to the system size of the linear-chain state. When f3 is large enough compared with 
the a-cluster size, the 1D-THSR wave function describes a dilute linear-chain gas where n a clusters move almost 
freely like a gas in the z direction. In the present calculation, the Ri integration is approximated by summation on 
mesh points in a finite box. The details are described later. 

The BB and 1D-THSR wave functions contain the total c.m.rn. In the present paper, the c.m.m. is not removed 
exactly because separation of the c.m. coordinate and intrinsic coordinates is technically difficult in calculation of the 
density matrix. In analysis of entropy in a system of la (one a) and that of an a cluster with a core, I keep the total 
c.m.m. as defined in the original form. For 2a, 3a, and 4a systems, I make a correction of the c.m.m. to eliminate a 
possible artifact from j3 dependence in the c.m.m. as explained later. 


B. Entanglement entropy 

I briefly describe entanglement entropy defined by the one-body density matrix. More detailed explanations are 
given in appendixes. 


1. Density matrix 


For a wave function | 4 h 4 )) f or an A-nucleon state, the one-body density matrix p^> is defined in the coordinate 
space as 


p (1 \ra\ r'a') = (T 1 A )|a’^(r , cr , )a(rer)|'I' (j4 )) 


(5) 
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where a'{ra) and a(ra ) are creation and annihilation operators of a nucleon at the position r with spin-isospin 
a = p t,P I, n t, n 4_. The one-body density matrix is regarded as the matrix element of the density operator p{^ for 
the wave function 

p (1) (ro-;rV) = {ra^^r'a'). (6) 

The diagonal element of the density matrix p^\ra) = p^\ra;ra) indicates the one-body density of a nucleons at 
r. Using the orthonormal single-particle bases {| Z)} that diagonalize the density matrix, the density operator and 
density matrix are written as 

4 1} = Ei *>pz (1) <ji. ( 7 ) 

l 

p (1) (rcr;rV) = ^2(ra\l}p[ 1) {l\r'a'}, 

l 

= ^Mr^Pi^tfir'v'), (8) 

l 

where (j)i{rcr) = (ra\l) is the wave function for the single-particle state |Z), and 

p\ 1] = (v^la^l^)}, (9) 

0 < p[ 1] < 1, (10) 

is the eigen value of the density matrix and indicates the occupation probability of the single-particle state |/}. The 
trace of the density matrix p U' equals to the particle number as 

A = Trp^ = ElE ~ E / P^ix^dr (11) 

l <T d 

P {1) (ra) = ^p^tfiraWiira), (12) 


2. Entanglement entropy 


The entanglement entropy is defined by the one-body density matrix as, 

= — Trp^ log// 1 ) = — EzE l°g Pi ^■ ( 13 ) 

i 


The entanglement entropy is zero if a wave function is a Slater determinant, because pj 1 ^ = 1 for occupied states 

and pi' 1 = 0 for unoccupied states. Thai is, the density operator p^ 1 satisfies {p^}} 2 = p ^ in the single-particle 
Hilbert space for a Slater determinant wave function [29(. 

In analogy to the expression D3J for the particle number by the a sum and r integral of the “local” density p ^ ( ra ), 
I define “local” entanglement entropy as follows, 


= E l s^\ra)dr, 
a d 

s (1 \ra) = E [-pf )lo gP ; (1) 
l 


(j)*(ra)(j)i(r(j). 


(14) 

(15) 


Here the factor [—p^logpj 1 )] is contribution of the single-particle state |Z) in S^\ and <j>*(ra)4>i(ra) means the 
density distribution in \l) and it is normalized as J2& f dr^>*(ra)<j)i(ra) = 1. Therefore, the local entanglement 
entropy s^(rcr) reflects spatial distributions of the important states 1 1) that contribute to the total entanglement 
entropy, whereas it is hardly affected by almost occupied states having p[^ ss 1. 
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3. Entanglement entropy for linear-chain a-cluster states 

For a system of independent Fermions, the wave function is given by a Slater determinant and it has zero entan¬ 
glement entropy, S'*' 1 '* = 0. In general, the entanglement entropy indicates how particles are entangled with other 
particles, and can be an indicator for many-body correlations. 

In a na-cluster state, the total wave function is spin and isospin symmetric and the one-body density is block 
diagonal with respect to er; p* 1 ) (rer; r'a')) = pA) (r;r')5 a(T >, where the reduced matrix p( 1 )(r;r / ) is independent to a. 
Therefore, I can discuss the density matrix and the entanglement entropy with the reduced matrix in the subspace, 


A = 4 n, 

(16) 

n = J p(r)dr , 

(17) 

p{r) = p^\ra, ra), 

(18) 

l— 1 

II 

» Co 

(19) 

S = J s(r)dr , 

(20) 

s(r) = s^^(rcr). 

(21) 


In the present paper, p(r), S, and s(r) indicate the density, the total entanglement entropy, and the local entanglement 
entropy, respectively, defined by the reduced density matrix. 

In the present ID cluster wave functions, single-particle wave functions for x and y coordinates are common for all 
nucleons and give no contribution to the entanglement entropy. Therefore, I discuss only 0 dependence of the local 
density and the local entanglement by integrating x and y coordinates as 


p{z) = J p( r )dxdy, 

(22) 

,(z) = S ,(r)ixdv ' 

(23) 


C. Model wave functions for na, 8 Be, and ”°Ne systems 

1. Wave functions for linear-chain 2a, 3a, and fa states 

For 1D-THSR wave functions of linear-chain 2 a, 3a, and 4a states, the R % integration is approximated by summation 
on mesh points at 1 fm intervals. To eliminate the /3 dependence of the c.m.m., I correct the c.m.m. by shifting the 

cluster position, Ri —»• R[ = Ri~ Rq of basis BB wave functions (4>g§) in 4>™g_ THSR (/3) with Rq = (Ri 4-b R n )/n. 

Consequently, the 1D-THSR wave function is approximated as 

{ " c>21 

~ f ■ ■ ■ ) R'n)i (24) 

i -1 P ) 

where mesh points of the summation are truncated in a finite box |i?'| < 12 fm. The correction of c.m.m. is equivalent 
to replacing the /3 dependent c.m.m. &g(/3) in the original 1D-THSR with 4 >g( 0) localized at the origin to eliminate 
the /3 dependence in the c.m.m., 

^1D-THSR = <W)*m t G9) ^ 4>&(0)<fr in t(/?), (25) 

sm cc exp j-jixS, - A.xl, - < 26 > 

where Xg is the total c.m. coordinate. 

2. Intrinsic wave function for s Be(0 + ) 

For 8 Be(0 + ), I use 2a-cluster wave functions. Funaki et al. have shown that 8 Be(0j*") is well described by the 
three dimension (3D) THSR wave function of 2a Si. In the present paper, I consider the spherical 3D-THSR wave 
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function, which can have 98% overlap with the exact solution of the 2a state for 8 Be(0^“) [24|. The definition of the 
3D-THSR wave functions is explained in Appendix [E] 

I use the P 2 -weighted 1D-THSR wave function of 2a, 

&L(&\P) = J dPP 2 exp{-|^j ^(ith = +R/2,R 2 = -R/2), (27) 

which is regarded as the intrinsic state of the spherical 3D-THSR wave function before the angular momentum 
projection as 

P J = 0 <b“(P 2 ;/3) « J dPexpj-^} = +^,R 2 = -|) = §^<& 2 g. THSR (/3), (28) 

where the factor 3> G (0)/i> G (/3) is the correction of the p dependent c.m.rn., and 4> 3 D _ T hsr(/ 3) is the spherical 3D- 
THSR wave function, which are described in Appendix [E] In the practical calculation, the R integration in Eq. 63 
is approximated by summation of mesh points at , R = 0,1,..., 24 fm. 


3. Intrinsic wave functions for 20 Ne(0 + ) and 20 Ne(1 ) 


For 2O Ne(0f) and 20 Ne(l 1 ), I use 16 0 + a-cluster wave functions. Zhou et al. have shown that these states of 
20 Ne are well described by 3D-THSR wave functions of 16 0 + a G. 7], In the present paper, I consider the spherical 
3D-THSR wave functions, which can have more than 98% and 99% overlaps with the exact solutions of the 16 0 + a 
states for 2O Ne(0 + ) and 20 Ne(l“) [24j |. 

In the present paper, I fix the center position of 16 0 at the origin by omitting the recoil of the 16 0 core and 
consider the P fe -weighted Gaussian distribution of an a cluster around the 16 O core. The intrinsic wave function of 
the 3D-THSR wave function of 16 0 + a is essentially given by the following separable form of the 1D-THSR wave 
function of a-(2a) on the z axis and the (0p x ) 4 (0p y ) 4 ho shell-model wave function as 

^ e (R k ;P) = A {$ho [(0 Px ) 4 (0 Py ) 4 ] $“ D ( _THSR(^ fc ; /3)} (29) 

^id-thsrJ P) = J dRiRl exp | -| $|“ B (Pi , P 2 = +e, P 3 = -e ), (30) 


where $ho [( 07 'a;) 4 ( 0 pj / ) 4 ] is the 8-nucleon wave function of the ( 0p x ) 4 (0p y ) 4 ho shell-model configuration at the 

origin, and ^iD-THSR^k’ ^ 1D-THSR a-(2a) wave function for an a cluster with the i? fc -weighted Gaussian 

distribution around the fixed 2a core at the origin, and e is taken to be an enough small value. ® 2 o Ne (R k ', P) is 
regarded as the intrinsic wave function of the spherical 3D-THSR wave function for 2O Ne(0 + ) and 20 Ne(l~) in the no 
recoil approximation as 


P J=0 $S t Ne (P 2 ; P) « y<iJl 1 exp|-M|$“o-a(Ji' = o,Jl 1 ), 

P^ =1 ^ e (R 3 ;P) « ydfli/2iyi M (Ai)exp|-M|$“o-a (J j/ 


$ 


b O-a , 
BB 1 


Ii'.Ih ; 



a 

Hi 


0, Hi), 


(31) 

(32) 

(33) 


Here is the ho shell-model wave function of the p-shell closure around the position R!. 

In the intrinsic wave function $ 2 o^ e (P fc ; /3), 8 nucleons in 0 p x and 0 p y orbits give no contribution to the entanglement 
entropy, and therefore, I analyze the 1D-THSR wave function of a- (2a), T^j^gR( R k ; 8 ), for 20 Ne. In the practical 
calculation, the Pi integration in Eq. (l30l) is approximated by summation of mesh points, Y^r-o ±i ±12 fm> an< ^ 
e = 0.02 fm is used. 


III. RESULTS 


I analyze the system size (/3) dependence of the entanglement entropy in ID cluster states of na and 16 0+Q! systems. 
Based on the analysis of the entanglement entropy, I discuss the delocalization of clusters in the linear-chain 3a and 4a 
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states, the 2 a state for 8 Be, and the 16 0+a states for 20 Ne, whose intrinsic wave functions are approximately described 
by the 1D-THSR wave functions, 4>"g_ THSR (/l), < ^ > 8Be(^ 2 ’^)’ an d ^iD-THSR(-^ fc ’ ^)> respectively, with optimum /? 
values. The optimum /3 values for the linear-chain 3a and 4a states are taken from Ref. [gj, and those for 8 Be(0^"), 
20 Ne(0 + ), and 20 Ne(l“) are reduced from the results for the spherical 3D-THSR wave functions in Ref. [24 ]. 


A. Analysis of one cluster 

I analyze entanglement entropy of a la system expressed by the 1D-THSR wave function < h[g_ THSR (/3), where an 
a cluster moves in the Gaussian distribution with the range /3 around the origin. The c.m.m. is not corrected. This 
system is regarded as an a cluster trapped in a ho external potential. In 4 >[q_ thsr (/3), a is the composite particle 
of 4 nucleons, but it is easy to mathematically extend the 1D-THSR wave function for a general composite particle 
consisting of Nf constituent particles ( Nf is the particle number in the composite particle), 

N f 

4% = II ^RlXcr, (34) 

(7 = 1 

where er is the label for Nf species of particles, which are not identical to each other. Figure[T](a) shows /3 dependence 
of the entanglement entropy S defined by the reduced density matrix for a a particle. In the small system size (/3) 
limit, S goes to zero because the cluster is localized around the origin in the 1D-THSR wave function. With increase 
of the system size /3, the entanglement entropy S increases. For a fixed /3, S is saturated with increase of the number 
of constituent particles (Nf). 



0 2 4 6 8 

P (fm) 


FIG. 1: (a) /3 dependence of entanglement entropy S for a cluster consisting of Nf particles for Nf = 2, 4, 8, and 16. S defined 
by the reduced density matrix for a fixed species a is shown, (b) Same as (a) but /3 dependence of exp(5'). 

If the one-body density is fragmented equally into to states, S = log to (to is the number of states). For instance, in 
the case that the one-body density is fragmented equally into two states, S = log 2 = 0.693. It is naively expected that 
the number of independent states that contribute to the entanglement entropy S is proportional to /3/b (the system 
size divided by the cluster size). Provided that the one-body density is fragmented equally into these states, exp(S') 
should have a linear dependence on the system size (3. As shown in Fig. [TJb), exp(5') is a almost linear function of fi, 
supporting the naive expectation. 

Spatial distributions of local entanglement entropy s(^) and density p(z) in the 1D-THSR wave function for the 
la system is shown in Fig. [2] In case of the system size as small as /3 = 1 fm, where an a cluster is well localized 
around the origin as shown by the density almost the same as that for a fixed a cluster, s(i:) is quite small. As /3 
increases, delocalization of cluster develops and s(z) increases in particular in low-density regions at the surface. s(z) 
is relatively suppressed at the high-density region near the origin. 






The analysis of one-cluster systems suggests that the entanglement entropy is enhanced in a low-density cluster gas 
having a relatively larger system size than the cluster size. The entanglement entropy S and the local entropy 5 ( 2 ) 
can be a good measure for the delocalization of cluster. 
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FIG. 2: Spatial distributions of local entanglement entropy s(z) and density p(z) in the 1D-THSR wave function for the la 
system with /3 = 1, 2, 3, and 5 fm. The density po(z) of a localized a for the /3 = 0 limit is also shown in the left top panel. 


B. Linear-chain states of na systems 

I analyze entanglement entropy of the 1D-THSR wave functions for the linear-chain 2a, 3a, and 4a states. Figure 
[3] shows the system size (/3) dependence of S. Similarly to the la system discussed previously, S increases as the 
system size /3 increases. In Fig. |4j entanglement entropy per a cluster, S/n, is plotted as a function of /?. If a- 
cluster motion is not affected by other a clusters, S should be proportional to the number (n) of a clusters. S/n is 
somewhat suppressed as n increases because of the Pauli blocking effect between a clusters, which comes from the 
antisymmetrization effect between nucleons in different a clusters. That is, delocalization of a clusters is suppressed 
because the effective system size for the a-cluster motion is reduced by the Pauli blocking effect. 

The Pauli blocking effect should be relatively strong in the inner high-density region rather than at the low-density 
surface region. To see the contribution of the low-density tail part of the a-cluster Gaussian distribution to the 
entanglement entropy, I calculate S of the “tail-cut” 1D-THSR wave functions with no outer tail, where the basis 
wave functions are truncated as |I?'| < /? by cutting off the basis wave functions in the |i?(| > /3 region for the tail 
component. S of the tail-cut 1D-THSR wave functions indicates the contribution of the inner part in the total entropy, 
whereas, the difference of S between the tail-cut wave functions and the original ones approximately corresponds to 
the contribution of the tail part. S of the tail-cut 1D-THSR wave functions shown by dashed lines in Fig. [3] is much 
smaller than that of the original 1D-THSR wave functions in 3a and 4a systems, indicating the significant contribution 
of the tail component in S. The tail contribution in the total entanglement entropy is larger than 50% in the /3 < 4 
fm region for 3a and in the /3 < 6 fm region for 4a. For the same /3 value, the contribution of the inner part in S is 
relatively small in 4a than in 3a because of the Pauli blocking effect. 

Figure [5] shows local entanglement entropy s(z) of 2a, 3a, and 4a for the system size /? = 2, 5, and 8 fm. In 
case of the small system size as /3 = 2 fm, s(z) is almost zero, which indicates no delocalization of a clusters. As 
/3 increases, s(z) increases and shows a broader spatial distribution than the density distribution. s(z) is relatively 
suppressed in the inner high-density region and enhanced in the low-density region at the surface. It means that the 
significant contribution in the entanglement entropy originates in more broadly distributed orbits than the system 
size. Comparing the result for a given /3 between 2a, 3a, and 4a systems, the larger n (the number of a clusters) 
system shows more suppression of s(z ) in the inner region than the smaller n system because of the Pauli blocking 
effect. 

Let us consider correspondence of the present result with the linear-chain 3a and 4a states in 12 C and 16 O predicted 
by Suhara et al. In Ref. [I], they have applied the generator coordinate method (GCM) to the linear-chain 3a and 4a 
states using effective nuclear forces to exactly solve ID dynamics of 3a and 4a, and obtained the energy minimum wave 
functions in the model space of linear configurations. It has been shown that the wave functions of the linear 3a and 
4a states obtained with the GCM calculations have large overlap with the 1D-THSR wave functions. The optimum 
























9 


2 

1.5 
cn 1 
0.5 
0 

0 2 4 6 8 

P (fm) 

2 

1.5 
m l 

0.5 
0 

0 2 4 6 8 

(5 (fm) 

2 

1.5 

m 1 

0.5 
0 

0 2 4 6 8 

P (fm) 





FIG. 3: j3 dependence of entanglement entropy S of the 1D-THSR wave functions for 2a, 3a, and 4a (solid lines). S of the 
tail-cut 1D-THSR wave functions is also shown (dashed lines). 



FIG. 4: Entanglement entropy per a cluster, S/n, plotted as a function of f3 of the 1D-THSR wave functions for 2a, 3a, and 
4a. 

parameters of the 1D-THSR wave functions are /3 = 5.1 fm and /? = 8.2 fm for the 0 + linear-chain states of 3a and 
4a, respectively. Looking at the result of corresponding wave functions, <f>3g -thsr(/ 3 — 5 fm) and 4>fg_ THSR (/3 — 8 
fm), shown in Figs. [3] and [5j it is found that the entanglement entropy is generated by the delocalization of a clusters. 
In particular, in the linear-chain 4a states with the large system size as /3 = 8 fm, s(z) is broadly distributed with 
significant amplitude. In both cases of the linear-chain 3a and 4a states, the tail contribution is significantly large 
as shown in the comparison of S between the tail-cut 1D-THSR and the original 1D-THSR wave functions in Fig. [3] 
In the 3a state with /? = 5 fm, the tail contribution is as large as ~ 50% of the total entropy. In Fig. [Gl spatial 
distributions of s(z) and p(z) in the tail-cut 1D-THSR wave function of 3a with /? = 5 fm and those of 4a with 
/3 = 8 fm are shown by solid and dashed lines compared with those in the original 1D-THSR wave functions shown by 
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dash-dotted and dotted lines. For the 3a state with /3 = 5 fm, the tail-cut 1D-THSR wave function shows the overall 
reduction of s(z) compared with the original 1D-THSR, indicating that the delocalization of a clusters occurs mainly 
in the low-density region at the surface. In other words, about half of the total entanglement entropy is generated 
due to the small difference in the tail component of the a distribution of the original 1D-THSR wave function from 
the tail-cut 1D-THSR one. In the 4a state with /3 = 8 fm, the tail-cut 1D-THSR wave function shows some reduction 
of s(z) at the surface region, but a significant amplitude of s{z) still remains in the inner region even after the tail 
cut. It means that, in the 4a state, the delocalization of a clusters occurs also in the inner region as well as at the 
surface region. 

I also demonstrate s(z) and p(z) in the localized cluster limit of the linear-chain 3a and 4a states given by the 
BB wave functions, where each a cluster is localized at a certain position Ri. They corresponds to the conventional 
linear-chain states. I set the a-cluster positions Ri on the z axis at equal intervals d = 3.8 fm in the 3a system and 
d = 4 fm in the 4a system so as to give density ( p(z )) peak positions similar to those of the corresponding 1D-THSR 
wave functions, -thsr(/ 3 = 5 fm ) and $ id-thsr(/ 3 = 8 fm )- In Fig. ED s(z) and p(z) in these BB wave functions for 
the conventional linear-chain states are compared with those of the 1D-THSR wave functions. s(z ) in the BB wave 
functions completely vanishes because S is trivially zero for a single Slater determinant. This result indicates that 
local entanglement entropy s(z) shows the prominent difference between the localized and delocalized wave functions 
for the linear-chain 3a and 4a states even though the difference in density is not so remarkable, in particular, in the 
3a state. It is concluded that both S and s(z) are sensitive to the delocalization of a clusters. 
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FIG. 5: Spatial distributions of local entanglement entropy s(z) and density p(z) in the 1D-THSR wave functions of 2a, 3a, 
and 4a systems with /3 = 2, 5, and 8 fm. 


C. a+a for 8 Be 

As pointed out by Funaki et a/., the exact 2a wave function for 8 Be(0j~) is well described by the 3D-THSR wave 
function Si- In the present paper, I use the I? 2 -weighted 1D-THSR wave function, c f>8 1 g e (/3), in Eq. (1271) as the 
intrinsic wave function of the spherical 3D-THSR wave function. The parameter /3 = 3.3 fm is reduced from the 
optimum parameter B = 4.77 fm with b = 1.36 fm of the spherical 3D-THSR wave function for 8 Be(0^) taken from 
Ref. [HI, using the relation B 2 = b 2 + 2 /3 2 and the scaling b = 1.36 fm —» 1.376 fm. (The parameter B here is 
originally labeled by “cr” in Ref. j24j].l 

Figure [7(a) shows /3 dependence of the entanglement entropy S of 4>8g e (/3). For comparison, S of the i?°-weighted 
1D-THSR wave function is also shown. At the optimized parameter /3 = 3.3 fm for 8 Be(0^), S = 0.6 is generated 
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FIG. 6: Spatial distributions of local entanglement entropy s(s) and density p(z) in (a) the tail-cut 1D-THSR wave function 
of 3a with /3 = 5 fm, (b) the tail-cut 1D-THSR wave function of 4a with /3 = 8 fm, (c) the BB wave function of 3a with the 
interval d = 3.8 fm, and (d) the BB wave function of 4a with d = 4 fm. s(z) and p(z) in the original 1D-THSR wave functions 
of 3a with (3 = 5 fm are also shown by dash-dotted and dotted lines, respectively, in (a) and (c), and those of 4a with (3 = 8 
fm are shown in (b) and (d). 


by the delocalization of a clusters. Spatial distributions of local entanglement entropy (s(z)) and density (p(z)) in 
$8B e (^ = 3-3 fm) are shown in Fig. [7(b). 
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FIG. 7: (a) /3 dependence of entanglement entropy S of the R 2 -weighted 1D-THSR 4>8^ e (/9) of 2a. That of the normal (R°- 
weighted) 1D-THSR of 2a is also shown, (b) Spatial distributions of local entanglement entropy s(z) and density p(z) in 
$^ e (, 8) with d = 3.3 fm for 8 Be(0+). 


D. le O+a for 20 Ne 

As shown by Zhou et al., 20 Ne(0^) and 20 Ne(l)") are well described by the 3D-THSR wave functions @, 0 ]. As 
approximated intrinsic wave functions of the spherical 3D-THSR wave functions for 20 Ne(0^) and 20 Ne(lj”), I use the 

R 2 - and f? 3 -weighted 1D-THSR wave functions, ^id.^hsrX-^’ 0) m Eq. (131)1) . respectively. For comparison, I also use 
localized cluster wave functions given by following parity projected BB wave functions with the fixed a-(2a) distance 
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d, 


Ki ah± ( d ) = $ 3 B c h(Ri=d,R 2 =+e,R 3 = -e) 

±<l>im(Ri = ~d, R -2 = +e, R 3 = -e). (35) 

<l ) gg 2Q ^ ± (c?) is given by the linear combination of two Slater determinants, and in the large d limit, it has S = log2. 
For the 1D-THSR a-(2a) wave functions, the optimum parameters /3 = 2.0 fm and /? = 2.6 fm for 2O Ne(0) 1 ") and 
20 Ne(lj~) are reduced from B = 2.39 fm and B = 2.97 fm, respectively for the spherical 3D-THSR wave functions 
with b = 1.46 fm in Ref. f24j using the relation B 2 = /3 2 + b 2 /2 and the scaling b = 1.46 fm —> 1.376 fm. For the 
positive- and negati ve-p arity projected BB wave functions, I use the optimum parameters d = 3.05 fm and d = 3.85 
fm taken from Ref. j24j, respectively. 

Fig. [5] shows j3 dependence of S of the 1D-THSR wave function, $iD 2 THSR(^ fe ’ /®) for 20 Ne, and also shows d 
dependence of S for the BB wave function, $gg 2 “' ) ’ ± (d) for the localized cluster wave function. For comparison, 
I also show the result of the i?°-weighted 1D-THSR wave function of an a cluster with and without the 2a core. 
In the i? 2 -weighted 1D-THSR wave function of a-(2a), S is zero at j3 —> 0 and rapidly increases in j3 < 1, and it 
gradually increases in the /? > 1 region. In the i? 3 -weighted 1D-THSR wave function, S is finite even at /3 —> 0 because 
^id-thsr^ 3 ! & ~t 0) is equivalent to the shell model limit wave function having 3 nucleons in the sd shell and a 
nucleon in the pf shell, which has S = — | log | — j log | = 0.562. With increase of /3, S of ^d-thsr(-^ 3 > @) becomes 

close to that of $id_thsr(-^ 2 > /?)■ I n case of the positive-parity projected BB wave function, S rapidly increases with 
the increase of d and it becomes constant S = log 2 in the d > 3 fm region. The negative-parity projected BB wave 
function shows small d dependence of S as S = 0.562 at d —> 0 and S = log 2 in the large d region. S = log 2 is 
generated in both the positive- and negative-parity projected BB wave functions with a large d 1 because the a-cluster 
wave functions are separated into two parts in z > 0 and z < 0, which have almost no overlap with each other. This 
indicates that even the localized cluster wave functions can have the finite entanglement entropy S = log 2 by the 
parity projection. In comparison of S between the 1D-THSR and BB wave functions, it is found that the entanglement 
entropy S ~ log 2 generated in the 1D-THSR wave functions of a-(2a) with j3 < 1.5 fm does not originate in the 
delocalization of the a cluster but it is understood by the effect of parity projection. With further increase of j3 in 
the 1D-THSR wave functions, S increases due to the delocalization of the a cluster and it becomes larger than log 2. 

Let us consider correspondence of the present result with 2O Ne(0f) and 20 Ne(lj"). As mentioned previously, 
^ > id- 2 thsr(^ 2 > P = 2-0 f m ) and $id_thsr(-^ 3 ) @ = 2-6 fm) correspond to the intrinsic wave functions of 2O Ne(0j 1 ') 
and 20 Ne(lj") described by the spherical 3D-THSR wave functions, respectively. $gg 2a ^ + (d = 3.05 fm) and 
Kf ah ~(d = 3.85 fm) correspond to the optimized BB wave functions for 2O Ne(0j h ) and 20 Ne(lj”). Figure [9] 
shows spatial distributions of local entanglement entropy (s(x)) and density {p(z)) in these wave functions. For 
2 °Ne(0j l "), $iD 2 THSR ( i?2 ;P = 2-0 fm) and $gg 2Q ^’ + (d = 3.05 fm) show similar distributions of s(z) and p(z) to 
each other. s( 2 ) almost vanishes in the —2 < z < 2 fm region because of the Pauli blocking from the core and 
it has amplitude only at surface regions. $ id-thsr (-^ 2 ’& = 2.0 fm) and $gg 2 “^ ,+ (d = 3.05 fm) have the total 
entanglement entropy S = 0.76 and S = 0.67, respectively. They are close to log 2, which is generated mainly 
just by the parity projection, indicating that the delocalization of the a cluster is weak in the intrinsic state of 
2 °Ne(0j 1 '). For 20 Ne(lj"), ^id.thsr ^ 3 ’@ = 2-6 fm) shows a broader distribution of s(z) in the outer tail part than 
that of < I ) bb 2q ^’ _ (^ = 3.85 fm). < I ) id_thsr (^ 3 )P = 2.6 fm) has the total entropy S = 0.90 larger than S = 0.69 of 
$ 3 g 2a -*’ (d = 3.85 fm). It means that, the entropy is somewhat generated by the delocalization of the a cluster in 
addition to S = log 2 from the parity projection. This additional entanglement entropy by the delocalization comes 
from a low-density outer tail of the wave function. That is, the delocalization of an a cluster occurs in the low-density 
tail part. 


IV. SUMMARY AND OUTLOOK 

In order to investigate delocalization of clusters in the ID cluster systems, I proposed a method of analysis using 
entanglement entropy defined by the one-body density matrix. I studied the entanglement entropy of the ID cluster 
states of na and 16 O + a systems, and discussed the delocalization of clusters in the intrinsic wave functions of the 
linear-chain 3a- and 4a-cluster states, 8 Be(0j 1 ’), and 20 Ne(Oj l ", lj”). I investigated the entanglement entropy of the 
1D-THSR wave functions and compared it with that of the BB cluster wave functions, and showed clear differences 
in the entanglement entropy between localized cluster wave functions and delocalized cluster wave functions. 
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FIG. 8: (a) /3 dependence of entanglement entropy S of the R fe -weighted 1D-THSR wave functions 4 > id_thsr(-R , s ; P) °f 0i-{2a) 
and the 1D-THSR wave function < 1?id-thsr(£0 °f 1° without the 2a core, and (b) d dependence of S of the parity projected 
BB wave functions < f > Bg 2Q ^’ ± (d). 






FIG. 9: Spatial distributions of local entanglement entropy s(a) and density p(z) in (a) < Fid_thsr(^' 2 > P) w hh P = 2-0 fm for 
2 O Ne(0j l "), (b) $ 3 g 2a ^’ + (d) with d = 3.05 fm, (c) 4>“q_ 2 ^s R (R 3 ; P) with P = 2-6 fm for 20 Ne(lj~), and (d) $bb 2 ^’ _ (c!) with 
d = 3.85 fm. 


I calculated the entanglement entropy of the 1D-THSR wave functions for linear-chain 2a, 3a, and 4a systems, and 
discuss the dependences on the system size (/3) and the number (n) of a clusters. With increase of the system size 
P, the entanglement entropy increases as the delocalization of a clusters develops. With increase of the number of a 
clusters, the entanglement entropy per a cluster decreases because of the Pauli blocking effect between clusters. In 
order to clarify the spatial regions where the entanglement entropy is generated by the delocalization, I defined the 
local entanglement entropy s(z). I found that the entanglement entropy is generated in the low-density part whereas 
it is relatively suppressed in the high-density part, indicating that the delocalization of clusters occurs dominantly 
in the low-density region but it is suppressed in the high-density region because of the Pauli blocking effect between 
clusters. 

Moreover, I discussed the delocalization of a clusters in the linear-chain 3a and 4a states predicted in 12 C and 16 0 
by Suhara et al. [s| based on the analysis of the entanglement entropy. In the linear-chain 3a state, the delocalization 
of clusters occurs dominantly in the low-density tail region while it is relatively suppressed in the inner region because 
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of the Pauli blocking effect. In the linear-chain 4a state having the significantly larger system size than the linear-chain 
3a state, the delocalization occurs in the whole system. 

I also analyzed the entanglement entropy of the l? fc -weighted 1D-THSR wave functions of 2a and a-(2a) systems, 
which correspond to the intrinsic wave functions of the cluster states in 8 Be(0^"), 2O Ne(0^), and 20 Ne(lj”). In 
8 Be(0f), the entanglement entropy is generated because of the delocalization of clusters. In 2O Ne(0^") and 20 Ne(lj~), 
the entanglement entropy is strongly suppressed in the inner region because of the Pauli blocking from the core. The 
entanglement entropy is generated at the surface region mainly because of the parity projection. In particular, in 
2 O Ne(0i"), the delocalization of the a cluster is weak and gives a minor contribution to the total entanglement entropy. 

The present result shows that the entanglement entropy is generated by the delocalization of clusters. The entangle¬ 
ment entropy is sensitive to the localization and delocalization of clusters and it can be a good measure to discuss the 
delocalization of clusters in ID cluster states. In the present paper, I investigated the entanglement entropy only in 
ID cluster wave functions, which are regarded as intrinsic states of 3D nuclear systems, and discuss the delocalization 
in the ID motion, i.e., the delocalization in radial motion of clusters. Zero entanglement entropy in the BB wave 
functions for localized cluster wave functions indicates that clusters are not delocalized but localized in the intrinsic 
systems. However, I should comment that, in realistic 3D wave functions, the ID cluster wave functions should be 
projected on to the angular-momentum eigen states. The angular momentum projection to 0 + states from the ID 
cluster wave functions usually causes the delocalization of clusters in angle motion. In that sense, the localized cluster 
wave functions may contain implicitly the delocalization in the angle motion even though the delocalization does not 
occur in the radial motion. The localization of clusters in intrinsic states, for which the delocalization occurs in angle 
motion by the angular momentum projection, is regarded as a kind of strong cluster correlations 

In principle, the present method with entanglement entropy can be extended to 3D cluster systems, but it is 
practically not easy to calculate entanglement entropy in 3D because high dimensional single-particle bases are needed 
in diagolization of the one-body density matrix in 3D. Moreover, as explained in Appendix[Fl the entanglement entropy 
in 3D systems should have strong dependence on the single-particle angular momentum j at the Fermi surface, because 
considerable entanglement entropy is generated in the angular momentum projection, which might wash out a pure 
contribution from the delocalization of clusters in the entanglement entropy. 
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Appendix A: Density matrix 

I describe density matrices. For details of the one-body density matrix, the readers are referred to, for example, 
Ref. [29| . An A-body density matrix for a wave function 4 >{ra 1 r 2 cr 2 , ■ • ■, tact a) for an A-nucleon system is defined in 
the coordinate space as 


P (A) (r 1 CT 1 ,r 2 cr 2 ,..., r A a A ; r\a' Xl r' 2 a' 2 ,.. ., r' A a' A ) 

= A!4-*(r' a[, r' 2 a ' 2 ,..., r' A (j' A )^(ncn, r 2 cr 2 ,..., r A a A ) 

= (r^a'^a){r' 2 a' 2 ) ■ ■ ■ a ] (r' A a' A )a ] (r A e a) ■ • • a t (r 2 cr 2 )a t (rK7i)|4'}, (Al) 

where a'(r<r) and a(rer) are creation and annihilation operators of a nucleon at the position r with the spin-isospin 
a = p i,n t, ^ i- 'P is normalized as ('F|'P) = 1. p^ is regarded as the matrix element of the A-body density 
operator p ^ = |\E , )('F|. The one-body density matrix is defined in the coordinate space as 

p^\ra\r'o') = ('F ( ' 4 ^|a t (r , cr , )a(rcr)|4'^ j4 ^}. (A2) 

pi 1 ) i s given by the trace of the A-body density matrix 

p (1) (ro-;rV) = A ^ / dr 2 • *. dr A ^*{r'a', r 2 (T 2 , ■ ■ ■, r A (T A )^{ra, r 2 CT 2 ,..., r A a A ) 

(T 2 ,...,0-A J 

= ATr 2 ---TrV A) , 


(A3) 
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which is also called reduced density matrix. The one-body density matrix is regarded as the matrix element of the 
one-body density operator p[y 1 for the wave function 


p( 1) (rcr; r'a') = {ra^W a'), (A4) 

and I get 

p { 1 \ra;r'a') = ^ ^ P {ra)p^ ^p* q {r'a) 

PQ 

= ^2(f(r\p)p { ^(q\r'a'), (A5) 

pq 

where 

Ppq = <^ (J4) |4 C PI^ (A) ) (A6) 

is the matrix element of the one-body density operator pY in arbitrary orthonormal bases. p^ 1 is a Hermitian 
single-particle operator and it has the form 


p* ] = ^2\p)ppq{q\- (A7) 

pq 

The density matrix can be diagonalized by a unitary transformation of single-particle bases 

(ny^DV = P^Sw, (AS) 

aj = ^2Dvic \,, (A9) 

l' 

where 

p\ 1} = <* (A) |aJa«|* (A) >, (A10) 

0 < p, (1) < 1 (All) 

is the eigen value of the density matrix and means the occupation number of the single-particle state |^) in the wave 
function In the coordinate space, the diagonal element p^^ra) = p^\ra,ra) of the density matrix is the 

one-body density of a = p f, p 4-, n t, n | nucleons at the position r, and it is expressed in the bases 1 1} as, 

P (1 \rcr) = ^2p\ i) (th(ra)(j)i{r(j). (A12) 

l 

The trace of the density matrix /b 1 ) equals to the particle number: 

A = Trp (1) = J2 Ppp = P^ = H / P W ( r °) dr - ( A13 ) 

p l cr 


Appendix B: Entanglement entropy 

The von Neumann entropy of an A-particle system is defined by the A-body density matrix as 

S (A) = -Trp (A ) logp (A) , 

and trivially 5^ A ^ = 0 for the pure state |\I/) because {p^Y = P^■ 

The entanglement entropy is defined by the one-body density matrix as 

aS'( 1 ) = — Tr/b 1 ) logp ( B = — pY log p\^. 

i 


(Bl) 


(B2) 
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Note that the entanglement entropy can be defined by the general A'-body density matrix for A! < A — 1 which is 
obtained by partial trace of the complete A-body density matrix for some degrees of freedom, but in the present paper, 
I only consider the entanglement entropy for the one-body density matrix. The entanglement entropy is zero if and 
only if a wave function |\!/^) can be written by a Slater determinant because p^ = 1 or 0 for a Slater determinant 
as pj ^ = 1 for occupied single-particle states and p ^ = 0 for unoccupied states. It is equivalent to the following 
theorem: A wave function is a Slater determinant if and only if the corresponding density operator satisfies 

{/3^} 2 = Py 1 in the single-particle Hilbert space |29]]. It means that, for the ideal case of an uncorrelated Fermion 
system that the wave function is given by a Slater determinant, the system has exactly zero entanglement entropy, 
S A) = o. S « is finite only if a system contains many-body correlations beyond the expression of a single Slater 
determinant. It indicates that the entanglement entropy can be an indicator for many-body correlations. 

In the present paper, I define local entanglement entropy and analyze spatial distribution of the entanglement 
entropy. In analogy to the expression of the particle number by the a sum and r integral of the density p^\ra) } I 
define the local entanglement entropy as follows, 


5 d) 


s^(rer) 


f s^^(rcr)dr, 

G 

EH'W 

I 


<t>i{rv)<t>l(r(r)- 


(B3) 

(B4) 


(f>i(ra)<j>i{ra) means the density distribution in the state 1 1) and normalized as f 4>* ( r &)4>i{ r <j)dr = 1, and the 

factor [— p|^ log is the contribution of the state |Z) in the total entanglement entropy S^ 1 ). Therefore, s^(ra) 
reflects spatial distributions of the important states |Z) that contribute to the total entanglement entropy. Note that 
s(z) is not quantity determined only by local information at the position z. 


Appendix C: Entanglement entropy for correlated and uncorrelated systems in a toy model 

Let us consider correlated and uncorrelated wave functions in a simple toy model of an A particle system. For 
simplicity, A particles are assumed to be distinguishable and stay on sites in a space. The number of available sites 
(single-particle states) is N s and I use the label k 3 for the jth single-particle states. Each particle can occupy one of 
single-particle states, kj (j = 1,..., N s ). A wave function 4/ of an A-body state is expressed by a linear combination 
of direct products as 

T(1,2 ,...,A) = ££-£ C*(q(i, Ot2-> • • • ? ^a)0q:i (l)0Q2 (^) * * * (A), (Cl) 

OL1 ot 2 (A-A 

where oti=ki,h 2 ,..., and C{a i,..., qa) is normalized as |(^(1, 2,..., A)|4'(l, 2,..., A))| 2 = 1. 

Let us first consider uncorrelated systems. If a state is an ideal state of independent particles, the wave function 
can be written by a simple product of single-particle wave functions 

$(1,2,...,A) = ^i(1)^2(2).-^(A), 

Ai’i) = £ Ci(a)4> a (i) 

cx=ki,k2,...,k]y s 

and I get = 0 because the one-body density operator for the ith particle is given as = |V>i)(V ! i| an d obviously 
satisfies {p^’*} 2 = (Here p^’ 1 is the reduced one-body density matrix defined in subspace for the ith particle, 

and the entanglement entropy is defined by p* 4 -*’®.) Let us consider the state of free particles in zero momentum 
that all particles move freely in the whole system with an equal weight. The wave function is given as 

1 A 

T(l, 2,..., A) = —(*) + ( t>k 2 {i) H-+ 4>k Ns (*)] • (C4) 

N s i —i 

The wave function has S^' 1 = 0. Another example is a “localized cluster” system of a cluster, where all particles are 
localized at one site kj to form a composite particle (a cluster) at kj. The wave function is given as 

A 

*(1,2,...,A) = If &,(*). 

2=1 


(C2) 

(C3) 


(C5) 




This localized cluster wave function also has zero entanglement entropy, S = 0. 
Let us next consider the following example of a strong correlation limit, 
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*(1,2,..., A) 


1 


n ^ w+n ^ (*)+••• n ^ 



(C6) 


where A particles form a composite particle, and the composite particle moves freely in the whole space. This is a 
strongly correlated system, where, if a particle is observed at a certain site, all other particles are always observed at 
the same site. This is a strong coupling limit of the spatial correlation and corresponds to a delocalized cluster wave 
function of a cluster. I can easily get the one-body density operator 


n b 






3 =1 


(C7) 


and the entanglement entropy S^’ 1 = \ogN s . 

Thus, the entanglement entropy indicates how a single particle is entangled with other particles. A localized 
composite particle has = 0 because there is no entanglement even though particles have some spatial correlation 
in a sense. In contrast, if the delocalization of a composite particle occurs, becomes finite and it is proportional 
to the logarithm of the number of sites (single-particle states) where the delocalization occurs. It indicates that the 
entanglement emerges because of the delocalization of the composite parcle. 


Appendix D: Calculation of density matrix for linear na states 


In the present paper, I calculate matrix elements of the one-body density operator for linear-chain a-cluster states 
by the expansion of localized Gaussian bases, 


4>R k = ftb 2 ) 3/4 exp 


-w^- R ^ 2 


(Dl) 


with Rk = (0, 0, Rk)- For simplicity, I choose a species of nucleons, for instance, p f and describe only the spatial 
part of a = p f nucleons because an na-cluster state is spin-isospin symmetric. I take Rk with 0.75 fm intervals as 
Rk = 0.75 j fm (j = 0, ±1,..., 15), and prepare an orthonormal basis set {</> p (r)} from the non-orthonornral Gaussian 
bases {4>R k (r)} (k = l,... , 35). For the bases | p), the one-body density matrix is written as 

Ppq = (M/^lctcpI*^), 

= J drdr'<f>*(r )p (1) (r; r')(/> q (r'), (D2) 

M r ) = ( r \P )• ( D3 ) 

By the diagonalization of p pq , I can get the diagonalizing bases |Z) by the unitary transformation of | p). 


Appendix E: 3D-THSR wave functions of 2a 

The deformed 3D-THSR wave function proposed by Funaki et al. [3:] is given as 


^3D-<ithsr(/^-L) Pz) — dRidR 2 exp 


-£• 

2 = 1,2 


R,. x + Ri 


ft 


C i > BB(-Rl' -^2 


(El) 
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The spherical 3D-THSR wave function of the case /3_l = / 3 Z = /3 is written as 


® 3 D-thsr(^) — / dR.1dR.2exp 

oc A exp < — 




i=l,2 

X? 


R? 


^(Ri,R 2 ). 


(E2) 


f 3 2 + b 2 /2 


oc A 


$ g (/ 3) exp l - 


X 2 


2/3 2 + b 2 




4>( a i )4>(a 2 ) 


$g(/9) oc exp < - 


2X g 

P 2 + b 2 /2 


(E3) 


where Xg = (.Xu + X 2 )/2 and X = X\ — X 2 , and $>g{P) is the c. m. motion. The wave function can be rewritten 
as 


®3 D-thsr(/3) oc In! [ dRexp 


$g(0) 


R 2 


R 


R 




(E4) 


Appendix F: Extension to 3D system 

The analysis with the entanglement entropy is applied to ID cluster systems in the present paper. In principle, it is 
able to calculate entanglement entropy also in 3D systems. However, to extract information of correlations from the 
entanglement entropy of 3D systems, one may encounter a problem from trivial correlation due to the total angular 
momentum projection. Let us consider two particles (not identical to each other) in a spatial orbital with the angular 
momentum l. The L = 0 state after the total angular momentum projection is given as 2 /TtI^ A 4 ) ® ^A 4 ) an< l it 

has the finite entanglement entropy S = log(2Z + 1). In the case that a cluster develops spatially, nucleons in a cluster 
have strong spatial correlations, which are generally characterized by the mixing of high l configurations, and therefore, 
S may reflect the many-body correlation in the cluster. However, one should take care that, even for a 0 + state with 
a single j 2 configuration in the jj coupling shell model the entanglement entropy is finite as S = log(2 j + 1) because 
of the trivial correlation by the total angular momentum projection. It means that the entanglement entropy strongly 
depends on j of the major shell and such a large contribution from the angular momentum projection could make 
it difficult to extract information of pure correlations beyond the jj coupling configuration from the entanglement 
entropy. 
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